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1. INTRODUCTION 

The fast progression of wireless communications network systems within the last few years has 
become a growing need to advance and efficient computerized tools for radio wave propagation prediction in 
various environments [1]. The massive usage of smart gadgets such as laptops, tablets, smart phones, and smart 
watches, has continuously contributed to the advancement of wireless communication technologies to 
accommodate the demand in the market. Therefore, a smart radio wave propagation prediction model plays an 
important role in the analysis of wireless communication systems especially for the emerging smart city design, 
upgrading the 4G mobile networks, internet of things (IoT), and the evaluation of 5G network system [2]. 
Radio propagation is a foundation to analyze and predict the signal details [3]. Direct measurement is one of 
the methods to evaluate the onsite performance, but it needs a considerable amount of time and very costly. On 
the other way, the analysis can be performed by a computerized simulation tool that offers more dynamic, low- 
cost, and more efficient [4]. 

In recent decades, radio frequency (RF) propagation models have been one of the major interesting 
research subjects due to the rapid development of new wireless communication system. Some propagation 
models have been proposed to predict the propagation paths of the radio signals and determining the optimal 
locations of the base stations [5], [6]. Every approach is determined by different parameters such as coverage 
area, bandwidth, and types of environments. For bandwidth, the models can be classified as either wideband 
or narrowband. The wideband models can be further divided into arrival times of rays, power delay profile or 
tapped delay line approaches [7]. At the same time, narrowband models can be grouped into different types 
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based on path loss and large-scale or small-scale fading. In term of surrounding conditions, there are two major 
classifications which is either outdoor or indoor propagation. For indoor environment, the propagation models 
are grouped into stochastic, empirical or site specific models. Empirical models are derived from certain 
propagation parameters such as the operating frequency, the distance between the transmitter (Tx) and receiver 
(Rx), and the height of the antenna [8]. These variables are assumed to follow the regular experimental data 
from measurement campaigns. They usually represent a set of equations derived from extensive field 
measurements, simple and efficient to use and suitable for environments with the same characteristics. The 
implementation of this model is easier as the major parameters are only the Tx-Rx distance and the operational 
frequency that contribute to a lower computational load. The main shortcoming of this model is lower accuracy 
as the propagation environment has been generalized and not well defined. For example, the empirical model 
for macrocells cannot be used for indoor picocells because of its range specific and not site specific. 

In the stochastic approach, radio signal is modeled as random variables [9]. Similar to the empirical 
model, the precision is also lower for signal prediction as the details of the propagation environment are not 
considered. On the other hand, for site-specific models, Maxwell’s equations are used in the radio wave 
propagation prediction. In these models, the precision of the simulation result can be improved. The simulation 
environment is well defined, especially in launching and tracing the ray, which makes the radio wave 
propagation to be more accurate and practical [10]. However, due to the increase in the detail of the 
environment features and ray propagation to be calculated, most of the site-specific models require higher 
computational load and simulation time, which contributes to the main shortcoming of this approach. In order 
to perform radio propagation analysis, the ray tracing (RT) which belongs to site-specific models is the best 
option to generate the practical result [11]. For indoor environment, ray tracing technique is the most capable 
technology to increase the accuracy of the simulation results [12]. 

The efficiency of ray tracing method depends on its capability to simulate the real surrounding of the 
tested scenario and trace the path traveled by each ray from the transmitting base station. However, the 
complexity of the simulation will become higher as more detailed items in the environment are to be considered 
to increase the accuracy. The computational load to trace the path and other propagation mechanisms will be 
increased significantly [13]. There are a lot of methods studied by many radio propagation researchers, but 
only a few concentrates on the analysis of the adaptive step size for ray launching that may contribute to the 
optimization in modelling the indoor radio wave propagation. In this study, our proposed adaptive three- 
dimensional (3D) ray tracing was analyzed together with renowned 3D ray tracing using the shooting and 
bouncing (SB) approach. All the modelings and simulations are applied using an in-house developed simulator. 
The actual measurement layouts and experimental parameters were used in simulations for result verification 
and validation. The operating frequency of 3.5 GHz was chosen for this study, which belongs to the sub-6 GHz 
band for 5G technologies [14]. For comparison, the 3D shooting and bouncing RT was selected which is the 
most used technique in radio propagation modeling [15], [16]. This technique has been further expanded and 
used extensively for modeling radio propagation in multipath channels [17], [18]. Based on our analysis of the 
simulation results, the adaptive 3D ray tracing methods produced an accurate propagation prediction that may 
benefit the design of the simulator to analyze radio frequency signal propagation. 


2. RAY CONCEPT 

In the ray concept, every ray follows the nature of sunlight that is assumed to propagate from a source 
point and extends infinitely along a straight line. The simple way to obtain the power level of radio signal 
propagation in free space is by using Friis’s formula as shown in (1). 


x = GeGr () () 


In (1), P}, P-, Ge, Gr, A, and r represent the transmitted power, power level, transmitter antenna gain, receiver 
antenna gain, wavelength, and separation distance between the transmitter and the receiver, respectively. 
However, in ray theory, Maxwell’s equations are applied extensively in radio propagation modelling. The 
electric field is represented by (2) [19]. 


E(F) = e(r) eio S0) (2) 


The symbol (F) denotes the magnitude vector, S (F) refers to propagation path, and the wavenumber, 
Po —~ at high frequency. The tracing of ray for radio propagation simulation can be implemented using certain 
principles. Most rays are assumed to travel in a straight path within the same medium before intercepting an 
object or obstacle [20]. The main obstacles to be considered for indoor conditions must be electrically larger 
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than one wavelength. The reflection coefficient for radio propagation is calculated by Fresnel equations using 
the wave impedance as (3), (4): 


Z2C0S0;—-Z1cosO¢ 


ro 3 

S Z2cos0i+Z1cos0t ( ) 
Z2C0SO¢—Z1 cos; 

lp = || P2008 MESerces ve (4) 
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From (3) and (4), both vertical and horizontal polarizations reflection coefficients are represented by I’, and 
Ip, respectively. The reflected and incident angles are denoted by 0, and 6;. The impedances of the first and 
second mediums are denoted by Z; and Z» respectively. The impedance value is related to the relative 
permittivity of the object, e; and the permeability ur for non-magnetic material of the building structures. 

A ray is modeled as a line-of-sight (LoS) path if it propagates directly from a source to a receiver. On 
the other hand, non-line-of-sight (NLoS) path is considered if a ray is hitting any object and reflected once or 
multiple times when propagating to the receiver. Equation (5) is used to compute the path loss (PL) for LoS 
and NLOoS radio propagation [19]. 


PL(f,d)[dB] = FSPL(f,1m) + 10nlogyo aa +X, (5) 


By referring to (5), n denotes the path loss exponent and X, refers to a zero mean Gaussian variable pertaining 
to the standard deviation ø. The free space path loss (FSPL) represents the path loss for a 1-m distance between 
Tx and Rx, and can be calculated by (6). 


FSPL(f,1 m)[dB] = 20nlogyy “4 (6) 


In this (6), f denotes operating frequency while c refers to the speed of light. If a ray is hitting a sharp 
edge, a cone of rays is generated due to the diffraction mechanism. The analysis of the diffraction mechanism 
is more complex. Here in this study, the diffracted rays were calculated based on the knife-edge approach in 
the ray tracing simulation. The power level of the diffracted ray is very low as compared to direct rays and 
normally considered only for the shadowed area in the simulation. By using the geometrical theory of 
diffraction (GTD), the Luebbers model can be used to implement diffraction in RT [21]. The diffracted power 
level Erp is calculated using (7). 


e skp’ | pr 
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The symbols Eo, Dip p', p, k, denotes the source point amplitude, diffraction coefficient, the 
distance between the diffraction point to Rx, the Tx-Rx separation distance, and the wave number respectively. 


l a andl’ zi are the vertical and 


horizontal polarization reflection coefficients. The symbol @ represents forward directions, while nrm — 
@ denotes backward directions and ¢' refers to reverse phase. Other symbols used in (8) refer to specific 
equations as follows: B = @ + ¢',a + (B) = 2cos2[(2nnN + —8)]/2], L = (pp'(pt+ p”), 


For non-conducting objects, (8) is used to obtain D, as below. The symbols I’ 
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3. PROPOSED ADAPTIVE 3D RT METHOD 
Adaptive 3D ray tracing is a method to simulate the propagation route of each ray using mathematical 
and propagation mechanisms and generates the valid paths from the transmitter to the receiver. Each valid path 
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comprises of some straight lines connected based on propagation mechanism such as direct path, reflection, 
and diffraction. The trajectory direction to launch a ray is given by the angle, measured between the horizontal 
plane and a reference point, and the sub-angle, which is the azimuthal angle on a vertical plane, measured from 
the same reference point. The objective of this adaptive RT is to identify rays from the transmitter to the 
receiver in a flexible manner but efficient. It consists of four main steps: pre-ray launch, adaptive step-size 
calculation, updated ray launching and determination of ray reception. 

— At first, the ray is emitted at a horizontal angle based on the initial step size configured initially in the 
simulator. When the first ray is launched from the Tx, the signal power is stored as reference power. The 
second ray is then launched at a vertical angle with the initial step size, which is also set in the initial 
configuration. 

— Secondly, a new step size is calculated based on a comparison between the power level of the current 
launched ray and the power level of the previous ray. The increment or decrement of the new step size is 
determined by the adaptive step size parameter which is set on the configuration page before the simulation. 
The calculation of the next step size is based on these three conditions; if the power level of the current 
signal is less than the previous signal, the next step size is increased by one adaptive step size. If the power 
level of the current signal is higher than the previous signal, the next step size is decreased by one adaptive 
step size. If the power level of the current signal is the same as the previous signal, the next step size remains 
the same. The increment of the step size is limited to the maximum allowable step size set in the 
configuration page. On the other hand, the decrement of the step size is also limited to the predefined 
minimum allowable step size. These two parameters are required to prevent the step size to become zero or 
a very big value. The same calculation will be applied for both horizontal and vertical angles during the 
launching of rays from the base station. 

— Thirdly, the next ray launching at updated angles is conducted after the increment by the updated step size. 

— Lastly, the rays are traced throughout the simulation using the propagation principles including diffraction, 
reflection, and transmission to identify the valid rays that are received by the Rx. 

The significant improvement of the proposed adaptive 3D RT technique is the introduction of adaptive 
step size which contributes to the flexibility feature on the ray launching angle. For the region with a high or 
significant power level, the resolution of the ray launching is higher by having a smaller step size. On the other 
hand, for the area with weak signal strength, the resolution of the ray is reduced by having a larger step size. 
Figure 1 demonstrates the process flow of the proposed adaptive 3D RT approach. The process flow is repeated 
until the whole area of interest is covered. 

The complexity of adaptive 3D RT is low, as more rays are launched to an area with higher signal 
power. Because of the smaller ray launching step size at this region, better accuracy is obtained with higher 
numbers of received rays with high power level. Moreover, less computational time is needed as very few rays 
are launched to region with a high probability of weak signals. In the existing technique, rays are emitted to 
every direction, which contribute to a higher computational load to perform the intersection tests and other 
propagation calculations for all rays without any consideration to reduce the number of insignificant rays with 
a very low power level at the Rx or lost rays. However, the adaptive approach is much better in reducing the 
computational load due to its adaptive nature of ray launching based on strength of signal power. In the 
simulation, ray propagation mechanisms are combined which include transmission, diffractions and reflections 
before the ray hits the receptor [22]. The signal power strength, Ep can be obtained by using (9). 


En =E (Qn) (Tio RinArin) (Zo Tin Atin) [=o DmnAamn)e 7” (9) 


The terms Ein (Qn), an, bn, Cn, Rins Tin, and Dmr are the initial power level, amount of reflections, amount of 
transmissions, amount of diffractions, reflection coefficient, transmission coefficient, and diffraction 
coefficient, respectively. The correlated spreading factors are denoted by Arin, Atjn) and Agmnn. Sn represents 
the traveled distance of a ray. Each Rx receives multiple rays due to direct path, diffraction, and reflection. At 
the receiver, total power level E total is the summation of power for each ray, which is given by (10). 


bn n —jkSn 
Etotal = Eno Ein (Qn) (TZ RinArin) (Tjo TrnAtjn)(Tlm=o DmnAdmn)e di (10) 


The total number of valid rays for each Rx is denoted by M. The transmitting signal power P,, and 
antenna pattern, Q, is used to compute Ein (Qn) as given in (11). 
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The expression Ey = ZP, G, refers to the power level at reference distance of 1-meter from Tx. The term 


No, Gt, Gm, and G m refers to intrinsic impedance, gain of the antenna, antenna directivity, and polarization in 
the direction of Qn, respectively. The p, is the distance between Q, and Tx. Vn is the voltage with respect to 
the Rx antenna and polarization which is derived using (12). 


A2GrnRo 


Al 


Vn = (En: rm) e? (12) 


The term 7, Grn, Ro, Õ rm, and e/® are the wavelength, Rx gain, Rx impedance, Rx polarization, and fixed- 
phase shift, respectively. The receiver power level is then computed using (13). The overall number of valid 
paths is represented by M. 
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Figure 1. Process flow of the adaptive 3D RT method 
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4. VALIDATION OF RAY RACING RESULT 

In this research, the ray tracing simulation results were validated with the actual measurement data. 
The simulation layout was constructed according to the actual building where the measurement data were 
collected. The ground floor of the wireless communication center (WCC) building in Universiti Teknologi 
Malaysia (UTM) at Johor Bahru Campus in Malaysia was chosen as the venue for the measurement campaign 
to gather experimental data. The external walls of the building were made of concretes and bricks. The interior 
partitions for the main rooms were made using 5 cm gypsum board and thick glass. The doors and windows 
were made of wood, plastic, and transparent glass. The dimension of the floor was around 21 m by 30 m. 

Figure 2 shows the Tx and Rx setup for the measurement campaign. For Tx, an Agilent E4438C model 
signal generator was used together with the transmitting horn antenna with horizontal polarization and a beam 
width of 49 degrees. The height of the antenna is 1.5 meters and the gain is 14.5 dB. For Rx, an Agilent E4440A 
model spectrum analyzer is used and connected with an omnidirectional antenna with the same height of 1.5 
meter. The receiving antenna gain is 7 dBi. The Tx was fixed at one location in the tested area. For each 
measurement, a continuous wave (CW) signal with frequency of 3.5 GHz with power level of 10 dBm was 
transmitted. The measurements of the received power at the Rx points were recorded by the spectrum analyzer 
starting from 1-meter distance between Tx and Rx. The Tx and Rx antenna heights were maintained at 1.5 
meters during the measurements, as an important factor in result accuracy [23]. 

For comparison purposes, the most used 3D conventional shooting and bouncing ray tracing (SBRT) 
technique and the proposed adaptive 3D RT method were selected in this simulation. The simulation results 
from both ray tracing techniques were analyzed and compared with actual measurement data. The layout for 
the RT simulations was constructed based on the actual WCC building layout using the in-house simulator. 
Figure 3 shows the visual representation of the simulation layout of WCC in 3D. The red dot is the Tx position 
and the blue dots are the tested Rx points. 


Figure 2. The Tx and Rx setup for measurement Figure 3. WCC layout in 3D for RT simulation 
campaign 


In the experiment, a sample of 20 Rx points which includes the LoS and NLoS scenarios was tested. 
The power level and path loss were recorded for each point. For both RT techniques, the reports on the 
simulation results were generated in the simulator database which consists of information on the power level 
for each Rx point and the path loss. 

Figure 4 shows the comparison graph of the power level at Rx points between simulation results and 
measurement data. The root mean square error (RMSE) was calculated using the standard formula to analyze 
the difference in power level between the measurement data and simulation results. The conventional shooting 
and bouncing ray tracing (SBRT) method showed a power level RMSE of 5.85, with respect to the experimental 
data. Among the 20 points, the receiver points Rx4, Rx6, Rx7, Rx9, and Rx20 expressed small differences in 
power level with the 3D SBRT approach. By using the adaptive 3D RT technique, as compared to measurement 
data, the power level RMSE was 2.55. For this method, most Rx points except Rx8 and Rx33 showed very 
small differences with the measurements. In any case, if the RMSE is small, the simulation result is very close 
to the practical implementation as shown in measurement data [24]. From the power level analysis, the 
proposed adaptive 3D RT technique showed a smaller RMSE value which is 3.3 dB lower than the existing 
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SBRT method with respect to the measurements. It was proven that the proposed adaptive 3D RT technique 
had generated a result that is closer to the measurements as compared with the conventional SBRT approach. 
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Figure 4. Power-level comparison graph between the measurements, 3D SBRT method, and proposed 
adaptive 3D RT method 


Figure 5 shows the comparison graph of the number of rays received using the proposed adaptive 3D 
RT and the conventional 3D SBRT approaches. As the data on this parameter were not available from the 
measurement, the comparison was analyzed between the adaptive 3D RT method and the 3D SBRT data only. 
From Figure 5, the adaptive 3D RT method showed greater number of successful rays detected at the Rx points 
as compared with the conventional 3D SBRT technique. In RT approach, a higher level of accuracy for indoor 
propagation modeling can be achieved if a larger number of rays can be detected at the receivers. It is justified 
by the smaller value in power level RMSE for the adaptive RT method. 
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Figure 5. Number of rays received by the receiver: comparison between the conventional 3D SBRT and 
adaptive 3D RT methods 


Figure 6 shows the comparison graph of simulation time for the conventional SBRT and adaptive 3D 
RT methods. In this case, both methods used the server with the same specification and configuration. From 
the findings, the conventional method required more computational time to perform the simulation compared 
to the adaptive 3D RT method. Since both methods were classified as site-specific propagation models, they 
can offer more precise radio signal predictions for a specific scenario than the stochastic or the empirical 
models. However, high computational time is identified as the main shortcoming for these 3D site-specific 
models [25]. In this finding, the adaptive 3D RT method saved about 44.4% of the time over the conventional 
method. It shows that the proposed adaptive 3D RT method is faster in generating the simulation results as 
compared with 3D SBRT approach, and thus mitigates the weakness of the site-specific models. 
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In this study, we developed an enhanced ray launching algorithm with adaptive and flexible 
characteristics. By implementing an adaptive step size approach to launch more rays in the region with higher 
possibilities of strong signals as compared with other lesser probability region, it has greatly reduced the 
computational load and time which become the key contribution to this method. By analyzing the results, the 
proposed adaptive 3D RT technique showed a better improvement on simulation time, power level and 
coverage in modeling the radio wave propagation for the indoor scenario. 
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Figure 6. Comparison on simulation time between 3D SBRT and adaptive 3D RT approaches 


5. CONCLUSION 

In this paper, we performed an indoor radio propagation modeling at 3.5 GHz by using the proposed 
adaptive 3D RT approach. This technique improves the ray launching algorithm by emitting more rays in an 
area with a high possibility of high power by adjusting the step size based on previously received power, which 
minimizes the computational complexity. The analysis on power level at multiple Rx points was compared and 
verified with measurement data. From the findings, the proposed adaptive 3D RT technique showed a smaller 
RMSE value than the existing 3D SBRT method with respect to the measurements. At the same time, the 
number of successful rays detected at the Rx points is higher by using the proposed method. In RT approach, 
a higher level of accuracy for indoor propagation modeling can be achieved if a larger number of rays can be 
detected at the receivers. This adaptive method also saved about 44.4% of the simulation time over the 
conventional method. It shows that the proposed adaptive 3D RT method is faster in generating the simulation 
results as compared to 3D SBRT approach, and thus mitigates the weakness of the site-specific models. 
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